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This paper provides a new proof of a theorem of Chandler- Wilde, 
Chonchaiya and Lindner that the spectra of a certain class of infinite, 
I— I random, tridiagonal matrices contain the unit disc almost surely. It also 

obtains an analogous result for a more general class of random matrices 
whose spectra contain a hole around the origin. The presence of the hole 
forces substantial changes to the analysis. 
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O 1 Introduction 

T — I 

I 

^ Over the last fifteen years there have been many studies of the spectral prop- 

erties of non-self- adjoint, random, tridiagonal matrices A, some of them cited 

^ in [HI [131 HH It has become clear that if all of the off-diagonal entries 

Aij with i — j = ±1 of the matrices concerned are positive, the almost sure 
Umit as — i- oo of the spectra of random N x N matrices subject to periodic 
boundary conditions can be quite different from the spectral behaviour of the 
corresponding infinite random matrix, [HI [THl HSl [16] . Indeed the limit in the 
first case can be the union of a small number of simple curves, while the second 
limit has a non-empty interior. 

Numerical calculations suggest that the situation is quite different if the off- 
diagonal entries have variable signs, but much less has been proved in this 
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situation, which is the one that we consider here. In a recent paper, [5], 
Chandler- Wilde, Chonchaiya and Lindner made important progress in deter- 
mining the almost sure spectrum of a remarkably interesting class of non-self- 
adjoint, random, tridiagonal matrices introduced by Feinberg and Zee in |13j . 
and sometimes called random hopping matrices, because the diagonal entries 
all vanish. Specifically they proved that, contrary to earlier conjectures, the 
infinite, tridiagonal matrix 



A. 
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has spectrum that contains the unit disc almost surely, [5]. The paper assumed 
that the entries c„ are independent and identically distributed with values in 
{±1}. 

In the present paper we assume that the entries c„ are independent and identi- 
cally distributed with values in {icr} for some fixed a G (0, 1]. We assume that 
the probability p that = a satisfies < p < 1; the corresponding probability 
measure on Vl„ = {±0"}^ is denoted by /x. The matrix Ac is identified with the 
bounded operator acting in the natural manner on £^(Z). 



In Lemma 26 we prove that 



Spec(Ac) C {A : 1 - a < |A| < 1 + a} 
by a perturbation argument. We also prove that 



Spec(Ac) C{x + iy: \x\ + \y\ < ^/2{1 + a^)} 

by obtaining a bound on the numerical range of Ac- There are currently no 
general techniques for identifying the precise forms of holes in the spectra 
of non-self- adjoint operators, and we have not done so here, but numerical 
calculations are consistent with the hypothesis that it is the intersection, H^j, 



of two elliptical regions as defined in (36); see the figures at the end of Section 



[7| Little is known about the part of the spectrum of Ac outside the unit disc 
even in the case a = 1, but numerical studies suggest that the boundary of the 
spectrum has a self-similar fractal structure in that case; O [T7] . 

The main result of [5], that the spectrum contains the unit disc almost surely, 
is for the case that a = 1, when there is no hole in the spectrum. It depends 



2 



upon the identification of a particular sequence c G f2i sucli tliat tlie equation 
Acf = Xf has a bounded solution / for every A G C such that |A| < 1. 

Our Theorem [7] rederives the main result of [5], in which cr = 1, but depends 
on a certain operator identity introduced in the next section. Our main re- 



sult, Theorem 27, that the spectrum of A^. contains that part of the unit disc 



which is not in Hfj, applies to all a G (0, 1). We give a second proof of this 



result in Theorem 34, by combining some results for a = 1 with bounds on the 
Lyapunov exponents of certain transfer matrices. Both proofs depend, addi- 
tionally, on results on the spectra of operators on £^(Z) which have different 
periodic structures on the positive and negative half-axes. They also both de- 
pend on explicit spectral calculations which we are able to carry out for certain 
operators Ac with c having arbitrarily large period. 

Our main results, as just stated, concern the spectrum of the (bi-)infinite ma- 
trix Ac- In a shorter final section we spell out implications for the spectra of 
the corresponding semi-infinite and finite matrices, illustrating these observa- 
tions with computations of the finite matrix spectra. In particular we show, 
by applying recent results of Lindner and Roch ^19j, that, unlike A^, the semi- 
infinite matrix has no hole in its spectrum for a G (0, 1), but contains the unit 
disc for all a G (0, 1]. 

Let £a denote the set of all c ^ that are pseudo-ergodic in the sense of [TO] . 
Precisely, c G £^o- if for every finite sequence b : {1, ...,n} — )■ {±0"} there exists 
m G Z such that br = Cm+r for all r G {1, ...,n}. Such sequences c are easy 
to construct without any reference to probability theory. The following facts, 
proved in [10], and rederived in [181 E] instance of the application of limit 
operator arguments, will be crucial in this paper. 

Proposition 1 Ifb,cESa then Spec{Ab) = Spec(y4c). Let denote this set, 
which is the main object of study in the paper. If c E Qa, then c E almost 
surely with respect to the measure /i. Finally 



Sa= [J Spec (A;, 



To describe a further result we establish, for G N and a G (0, 1] let ttjvo- 
denote the union of Spec(y4c) over all c G fio- that are periodic with period 
< A^. Let 

TToo.a = [J T<'N,a- (1) 

One obvious implication of the above proposition is that 

TTocCT C Sa- (2) 
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As is well-known, the set ttat^o- is the union of eigenvalues of N x N matrices. 
(Precisely, it is the union, over all sequences c and all |a| = 1, of the eigenvalues 
of the matrix Ai^"^*^^^ defined in (51) below; see (40) and |TT]. For another. 



equivalent characterisation see Lemma 20 ) This simple observation is useful 



in that it provides a method for computing what prove to be large subsets of 
Scr, and will be one component in our arguments. 

An interesting question is whether tTqco- is dense in 5*0-. We do not answer this 
question one way or the other, but our method of proof of Theorem [7| showing 
that the unit disc is a subset of 5*1, as a by-product, and with some additional 
argument, leads to a proof that 71^0,1 is dense in the unit disc (Theorem 15). 

For the sake of simplicity we will, throughout the rest of the paper, omit the 
subscript a in our notations if a = 1. We use N and Z+, respectively, as our 
notations for the sets of positive and non-negative integers. 



2 An abstract theorem 



In this section we present an abstract theorem that might be interesting in 
other contexts. It will be applied in Section |3| 

Let y4 be a bounded linear operator acting on the Hilbert space T-L and let 
Ti = Tie® 'Ho be an orthogonal decomposition of Ti. 

Lemma 2 // AiT-Lf.) ^ Ho and A{l-io) ^ 'He then 'He and 'Ho are invariant 
under the action of A^ . If B is the restriction of A^ to Tie and M is the 
restriction of to 'Ho then 

Spec(A2)\{0} = Spec(5)\{0} = Spec(M)\{0}. (3) 

// A is invertible then 

Spec{A^) = Spec(5) = Spec(M). (4) 



Proof The decomposition "H = 'He © 'Ho allows one to write the operator A in 
the form 

■ X 
Y 



A = 

where X : Ho — ^ and Y : "He ^ 'Ho- Therefore 
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This imphes that B = XY and M = YX. The second identity in ^ follows 
by some simple algebra that holds for any pair of bounded operators X and 
Y, and the first identity is a trivial consequence. 

If A is invertible then ([s]) implies that B and M are also invertible; therefore 
Q is equivalent to g. □ 

Theorem 3 Let l-i = i'^(Z), let Tie be the closed subspace of sequences whose 
supports are contained in the set of even integers, and let Tio be the closed 
subspace of sequences whose supports are contained in the set of odd integers. 
Let A be a bounded operator on l-i whose matrix satisfies Ar^s = for all r, s 
such that \r — s\ 7^ 1. Then AiTie) C "Ho O'lT'd AiTio) C "Hg. Moreover the 
identities 

Spec(A2) = Spec(S) = Spec(M) 
are valid in either of the following two cases. 

1. \Ar^s\ = 1 for all r, s such that \r — s\ = 1; 

2. There exist constants /3, 7 such that < /3 < 7 < 00 and |Ar,s| ^ P if 
r — s = 1 while |^r,s| > 7 if r — s = —1. 

Proof 

Case 1. An elementary calculation establishes that there exists a sequence 
/ : Z — 7- C such that Af = 0, |/2„| = 1 for all n and /2n+i = for all n, so 
that A and B are not invertible viewed as operators on £°°(Z), and thus not 
invertible as operators on £^(Z) (see e.g. Theorem 2.5.2]). So G Spec(y4) 
and G Spec (5). Similarly there exists a sequence / : Z — )■ C such that 
Af = 0, |/2n+i| = 1 for all n and f2n = for all n. Hence G Spec(M). The 
result follows by combining this with ([3]). 

Case 2. The operator Al associated with the matrix 



Ar^s if r — s = — 1, 
otherwise, 



is invertible and satisfies < 7 ^. The operator Ar = A — Ai satisfies 

\\Ar\\ < p. Therefore A is invertible with 



\A-'\\ = \\A-^\I + ArA-^')-'\\ < ^ ' 



1-/3/7 7-/3' 

The proof is completed by applying Q. □ 
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3 The case a = 1 



The following lemma was noted in [3]. 

Lemma 4 If c & Q then Spec(Ac) is invariant with respect to both of the maps 
A — )■ A and A—)-— A. If X E S then A and iX lie in S . Hence S is invariant 
under the dihedral symmetry group D2 generated by these two maps. 

Proof The invariance of Spec (Ac) under complex conjugation follows directly 
from the fact that Ac has real entries. If D is the diagonal matrix with entries 
D^^^ = {-iy for all r G Z then DA^D-^ = iA_c, so 

Spec(A,) =2Spec(A_,). (6) 

Iterating this identity yields Spec(y4c) = — Spec(y4c). This proves the first part 
of the lemma. The second part follows once one observes that c G £^ if and 
only if — c E S. □ 

The formulae in ([T]) are related to those in [5l Prop. 2.1], in a way that we will 
make explicit in Section |4} However, nothing resembling the following lemma 
appears in [S]. 

Lemma 5 Given b E Q, let c = r+(6) E Q be the unique sequence satisfying 

Co = 1, C2„ + C2„+l = 0, C2„C2„_1 = 6n, (7) 

for alln E Z. Then A^ is unitarily equivalent to A^QMb acting in ^^(Z)©^^(Z), 
where 

iMbf)n = -fn^l + (C2n+1 + C2n+2)/n + fn+1 (8) 

for all f G £^(Z). Moreover 

Spec{Al) = Spec(Afe) = Spec(Mb). 

Proof One may write {Acf)n = Cnfn-i + fn+i for all n G Z, or equivalently 
Ac = VcR + L where = = /„_i and {VJ)n = c^fn for all 

Therefore 

Al = VcRVcR + LVcR + VcRL + L'^ 
= XcR^ + Yc + L'^ 

where and are the diagonal matrices with diagonal entries 
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The operator has two invariant subspaces 



ne = {fe f{Z) : f2n+i = for all n E Z} 

and Tio = f'iZ) Q Tie- After an obvious relabeling of the subscripts, the 
restriction of A^ to "He equals At, while the restriction of A^ to Ho is equal to 
Mb, as defined in ([s]). The final statement of the lemma is now an application 
of Theorem [3} case 1. □ 

We will exploit extensively the formula Spec(A^) = Spec(A;,) which appears in 
the above lemma. The equation Spec(Ab) = Spec(Mb) will not play a role in 
our subsequent arguments, but makes an intriguing connection between spectra 
of rather different tridiagonal operators. Extending this connection slightly, for 
b E Q define c = r+(6) and by 

{Mbf)n = /„_1 +i"(c2„+l + C2n+2)fn + /n+1, 

for all / G £^(Z). Then, arguing as we do above to show ([g]), we see that 

Spec(Mb) = 2Spec(Mfc). 

In particular, in the case h E S when, by Lemma |4| zSpec(Mb) = i'^Y>^c[Ai) = 
Spec(y4b), we see that 

S = Spec(A) = Spec(Mb). 

Thus, in studying S, we are studying both the almost sure spectrum of the 
infinite hopping-sign matrix Ab with respect to the measure /i, and the almost 
sure spectrum, with respect to the same measure, of Mb, a discrete Schrodinger 
operator with a particular, complex random potential. 

In the next lemma we define the square root of any non-zero complex number 
to be the root whose argument lies in (— 7r/2, 7r/2]. 

Lemma 6 IfbEil and c = T^{b) then X E Spec{Ab) if and only if ±^/X both 
lie in Spec(y4c). If X E S then ±\/X both lie in S. 

Proof Lemma |4] and Lemma [5] imply that the following statements are equiv- 
alent. A E Spec(y4fe); A E Spec(yl^) = (Spec(74c))^; either or — lies in 
Spec(74c); ±V^ both lie in Spec(y4c). 

If A G S* and b E £ then A G Spec(yl5) by Proposition [l] Lemma [S] implies that 

A G Spec{Al) = (Spec(A,))^ C S^. 

Therefore either \/X or —VX lie in S. The proof is completed by applying 
Lemma HI □ 
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Theorem 7 The set S contains 

[j [0,2^/2"]. (9) 

neZ+, re{0,...,2"+2} 

Hence S contains the unit disc in C. 

Proof For n = the theorem states that 

[0,2]x{l,2,-l,-z}c5. 

This foUows by combining Lemma |4] with direct calculations of Spec(ylc) when 
On = 1 for all n E Z (in which case Spec(y4c) = [—2, 2]) and when c„ = — 1 for 
all n E Z (in which case Spec(y4c) = z[— 2,2]). For larger n the first statement 
of the theorem follows by applying Lemma[6]inductively. The second statement 
is now a consequence of the fact that the set (jol) is dense in the unit disc. □ 



4 The maps r± 

A crucial role has been played in the proofs above by the nonlinear map r+ on 
Q introduced in Lemma [5} and this map will be key to the arguments that we 
make throughout this paper. And in fact a sequence which is almost a fixed 
point of r+ (in a sense made precise below Lemma [s]) is central to the proof 
of Theorem [T] in [5] , though the proof is quite different and no mapping r+ 
appears in [5]. 

The relationship between the above proof of Theorem [T] and that in [S] is 
clarified to some extent by the following. Building on the definition of r+ 
made above, let us define maps T± : fl ^ fl hj T±{b) = c where 

Co = ±1, C2n + C2n+1 = 0, C2„C2n-l = ^n, (10) 

for all n G Z. We also define the space inversion symmetry 6 — )■ 6 by 6„ = 
for all n G Z. 

Lemma 8 IfT±{b) = c then T^{b) = c. In particular T±{c) = c if and only if 
r=p(c) = c". Each of the equations T±{c) = c has exactly one solution. 

Proof Let c = r+(6) and d = r_(6). Then d^ = —1, d2n + c?2n+i = and 
d2nd2n~i = bn = &i-n for all w G Z. Therefore do = di = 1. Also 

d2n+l + d2n = '^l-(2n+l) + '^l-2n = d-2n + '^l-2n = 
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and 

d2nd2n-l = dl-2ndl-{2n--l) = C^2(l-n)-lC^2(l-n.) = = &n 

for all n G Z. Therefore d = T^(h) = c and d = c. 

The proof that c = T_{b) implies d = r+(6) is similar. The other statements 
of the lemma follow immediately. □ 

This paper and [5] use three different special sequences. The sequences c± are 
defined by r-|-(c±) = c±. It follows directly from their definitions that c+,o = 1 
and = —1 while c_,o = ~1 and c__i = 1. However 

for all n 7^ 0, 1. The paper [5] uses the sequence Cg such that Ce,o = Ce,i = 1, 
while Ce,n = c±^n for all other n. Because of the space inversion symmetry the 
use of c+ or c_ in any proof is really a matter of convenience. 

We now turn to the solution of the equation AcU = Xu where m : Z — > C is an 
arbitrary sequence. The eigenvalue equation is equivalent to the second order 
recurrence equation 

Lemma 9 Suppose that c & Q and Cn = Ci_„ for all n G Z; that Un+i + 
CnUn~i = Xun foT sovfie A G C aud all n & 7a and Mq = 0, Mi = 1; and that 
Un+i + CnUn-i = Xun foT all u ^ Z and Uq = 0, Ui = 1. Then = \u^n\ for 
all n ^ Z. In particular is bounded as n ^ oo if and only if Un is bounded 
as n ^ —oo. 



Proof If one puts Vn = U-n then 

Cn+lVn+l + Vn-l = C_„M_„_i + = Am_„ = XVn (H) 

for all n G Z. Define a : Z — )■ {±1} by = 1 and an/an-i = Cn for all n G Z. 
If one now puts Wn = anVn for all n G Z then (11) implies 

Wn+l + CnWn-1 = XWn 

for all n G Z. Since wq = vq = uq = it follows that there exist 7 such that 
Wn = 7Cn for all n > 1. But |c„| = 1 for all n, so one obtains I7I = 1 by 
evaluating this identity for n = 1. Therefore \un\ = l^nl = \vn\ = \iJ--n\ for all 
n> 1. □ 



Corollary 10 Letu+, U-, Ue Z ^ C be the solutions of Un+i + CnUn-i = Xu 
for all n G Z subject to uq = and ui = 1, if c is put equal to c+, c_, c, 
respectively. Then U-^^n = Ue,n for all n E Z. Moreover = \ue,n\ for all 

neZ. 



n 
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Proof The first statement is proved by an elementary computation. For the 
second we use Ce = and c_ = c^. Lemma [o] now yields 

1^— ,n| 1'^+,— n| 1*^6, —n| ["^einl 

for all n G Z. □ 

The main step in the proof of Theorem [T] in |5] is contained in the follow- 
ing proposition (we quote here the parts of [SI Proposition 2.1] which we use 
immediately or later in section Isj) . 



Proposition 11 Let he defined as in Corollary 10 and define Pij E Z for 
by the formula 

i 

Ue,i = J2p^'^^^~^ 

with pij = if j > i. Let Y denote the set of {i,j) G such that pij ^ 0. 
Then Pij G {0, 1,-1} for all i, j and {i,j) E Y if and only if one of the 
following holds. 

(1) 1 = 3 = 1; 

(2) i and j are both even and (i/2,j/2) G Y ; 

(3) i and j are both odd and {{i + l)/2, (j + l)/2) G Y ; 

(4) i and j are both odd and {{i - l)/2, (j + l)/2) G Y. 

The following result is an immediate corollary of this proposition and Corollary 



10, which together imply that \ue^i\ < (1 — |A|) ^ for i G Z and |A| < 1. 



Theorem 12 As in Corollary\l(^ /ei : Z — > C be the solution of Un+i + 
CnUn-i = Xun for all u E Z subject to uq = and Ui = 1, with c = Ce- Then 
Ue E i°°{Z) for |A| < 1, so that Spec(y4c^) contains the unit disc. 



Since Corollary [10] has shown that m+ = Me, it is clear from Theorem 12 that 
Spec(74c^) also contains the unit disc. In fact this is precisely its spectrum. 

Theorem 13 If is the unique solution o/r_|_(c) = c then 

Spec(A,J = {z : \z\ < 1}. 

Proof It remains only to show that Spec(Ac^) C {z : \z\ < 1}. If A G 
Spec(y4c_^) then repeated apphcations of the first part of Lemma [6] yield A^" G 
Spec(74c^) for all n > 1. Since the spectrum is a bounded set, it follows that 
|A| < 1. □ 

We will (rather arbitrarily) focus on the mapping r+ rather than r_ in the 
remainder of the paper. The following lemma, which shows that the set of 
periodic sequences is invariant under the action of r_|_, will play a key role. 
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Lemma 14 If b ^ Q is periodic with period N, i.e. hn+N = bn, n & Z, then 
c = r_|_(6) is 4:N -periodic. Conversely, ifbEQ,c = r+(6), and c is 2N -periodic 
for some G N, then b is N -periodic. 

Proof First note that, if c = r+(6) and one defines c G ^2 by c„ = C2n, G Z, 
then 

C = r+(6) <^ (Co = 1, Cn = -bnCn-1, C2„+l = -C„, U G Z). (12) 

Therefore 

n 

Cm+n — Cm ( 1) ^m+j (13) 

i=i 

for all m G Z and n G N. If 6 is A^-periodic, then 

i=i i=i 
for all m G Z. Therefore c is 4A^-periodic. 

Conversely, if c = r+(6), for some 6 G fi, and c is 2A^-periodic for some G N, 
then c is A^-periodic and, from (12), it follows that b is A^-periodic. □ 

To illustrate the above lemma, define c^,c+ G ^2 by c~ = —1, c+ = 1, for 
n G Z, and define the sequences c*^™'+), c^"*--) g for m = 0, 1, by 

c(o,±) = c{m,±) ^ r+(c('"-i'±)), m G N. (14) 

Then explicit calculations of the action of r+ yield that c'-^'^-* = T^(c~^) is 
4-periodic (but not periodic with any smaller period), with cL^j^"* = Cq^'^"* = 1, 
c^^'"*"^ = C2^'^^ = —1. On the other hand, c*^^'~^ = r_|_(c~) is 2-periodic (and so 
also 4-periodic), with oil' •* = (—1)" for all n G Z. 

Both these calculations, of course, are consistent with the lemma, which implies 
that c*^"^'^^ is A^-periodic with A^ = 4*^, so that, using the notation ([l| (dropping 
a given that a = 1), 

Spec(74^(m,±)) C 7r4m, m = 0,l,.... (15) 

Although we do not have an explicit formula for the sequences c^'"'^^ it is easy 
to compute Spec{A^(m,±)). By Lemma |6l if c = r+(6), then 



Spec(AJ = {±v^ : A G Spec(A)}. (16) 
The proof of Theorem 17] begins with the observation that Spec(Ac+) = [—2,2] 



and Spec(Ac-) = i[— 2,2]. Combining this observation with (16) we easily 
prove by induction that 

Spec(A,(™,+)) = {re^'^/2'" . q < r < 2^2'", j e {0,...,2™+i - 1}} (17) 

11 



and 

Spec(A^(™,-)) = e''*/^'"^' Spec(A^(™,+)). 



Combining equations (15), (17) and (18), we see that we have shown that 



{r e-^^/^'" : < r < 2^/^-+^^ ^ ^ |o, 2^+^ - 1}} c 7r4^, m = 0, 1, .... 

Thus we have shown the following modification of Theorem [7] which, of course, 
by ^ , has Theorem [T] as a corollary. 

Theorem 15 The set -k^o contains the set and so is dense in the unit disc 
in C. 

We know Spec(y4^{m,±)) explicitly, but do not have explicit formulae for the 
sequences c^^'^\ However we can show that c''™'^^ converges pointwise to the 
sequence c+, the unique fixed point of r+, as m — )• oo. This is the content of 
the next two lemmas. We omit a proof of the first of these lemmas which is an 
easy consequence, by simple induction arguments, of the definition of r+. 

Lemma 16 Ifb&Q and c = r^{b), then cq = c+^o O'lT'd C\ = c+,i. //, for some 
N G N , bm = c+^m for m = 1, ...,N , then also Cm = c+,m for m = 2, 2N + 1. 
//, for some N e Z+, b_rn = c+ for m = 0,1, ...,N, then 6.^ = c+ _m for 
m = 1,2,...,2N + 2. 



Lemma 17 Let 6 G fi, and define c*^"^ G for n E ISi by c^^^ = r+(6) and 
^{n+i) ^ r+(c(")), n G N. Then, for n G N, 

c^^) = c+,m, m = 2-2'', 3-2^, 2" - 1, 

so that c*^"^ — )■ c+ pointwise and A^{„) converges strongly to Ac^ as n — )■ oo. 
Further, 

Spec(A,(„)) C {A : |A| < 2i/2"}. 

Proof The first equation follows by induction from Lemma [161 The second 



equation follows by induction from (16) and the trivial bound that Spec(74;,) C 



{A : |A| < 2}, which holds for all b e Vi. □ 



5 The mapping T^j^^ 

For the rest of the paper we consider operators for which the coefficients 
c„ take values in {±cr}, where < a < 1; that is, in the notation we have 
introduced in the introduction, we assume that c G Q^, for some a G (0, 1]. 
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The mapping that we have introduced continues to play an important role. 
We extend the mapping so that it operates on fi^-a, defining, for cr e (0,1], 
r^_+ : n^2 by 

r.,+ (c) =ar+((T-2c). (19) 
In other words, for b G c = T(j^+{b) is the unique sequence in Q„ satisfying 

Co = (T, C2n + C2n+1 = 0, C2nC2n-l = K- (20) 

Main properties of the mapping Fg- ^ for our purposes are contained in the 
following extension of Lemma |6j We will need to refer to a number of circular 
annuli, and use in this lemma and subsequently the notation 

|a,6l = {A : a < |A| < 6}. (21) 

Lemma 18 If h E VL„2 and c = Tcr^+ib) G Qa, then 

(A G Spec(A)) ^ (±v^G Spec(^)). (22) 

Hence 

iXeS^2)^{ i^A G ) (23) 

and 

{la,bjCS^2)^ila'/',b'/^CS,). (24) 

Proof We modify the calculations in Section [3j Lemma |4] is valid as it stands. 
In Lemma [5] we assume that b G il„2, and define c G by c = T^^^(b), and 
apply Case 2 of Theorem |3] in place of Case 1. This leads to the conclusion 



Spec(y4^) = Spec(y4b) as in Lemma pj (23) and (24) follow by choosing b G £"0-2 



and using Proposition [TJ □ 



6 Periodic and paired periodic operators 

To prove our main theorem we need results on operators Ac on £^(Z) that 
have one periodic structure for n > and another for n < (which we term 
paired periodic operators). The essential spectrum of such an operator is the 
union of the essential spectra of the periodic operators involved, which may be 
calculated explicitly using their Bloch decompositions. 

There may also be substantial inessential spectrum, in particular, open subsets 
of the spectrum where Ac — XI is Fredholm but has non-zero index. These parts 
of the spectrum (and the corresponding values of the index) can be computed 
by application of general results for block Toeplitz operators, which have been 
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developed to a high degree of sophistication; see [H 131 E], H] and the references 
therein. We need only a small part of this theory, and it is easy to develop this 
from first principles. We do this in a short Lemma [25] below, inspired by earlier 
analysis in [121 El HO], and particularly [101 Theorem 12]. Both the proof of 



Lemma |25| and the effective application of this lemma to prove Theorem [27 
depend on the next two lemmas which describe properties of the spectra and 
eigenfunctions of periodic operators. 

We assume throughout this section that the parameter a G (0, 1). 
Lemma 19 Lei 

, , Re(r)2 Im(r)2 ^ ^ 

•where r G C and — 1 < 7 < 1. Then the quadratic equation 

z'^ -Tz + -f = (26) 

has a solution satisfying \z\ = 1 if and only z/ $ = 1. // $ < 1 then both 
solutions satisfy < 1. // $ > 1 then one solution satisfies < 1 and the 
other satisfies \z\ > 1. 



Proof For G R, -2 = e*^ is a solution of ([26[) if and only if 



cos(fc') = , sm^fc') — 



1+7 1 - 7 

so that ( [26] ) has a solution satisfying \z\ = 1 if and only if $(r, 7) = 1. 
The set [/ = {r G C : $ < 1} is connected and contains the origin. Since the 



solutions of (26) depend continuously on r, and both solutions satisfy \z\ < 1 
if r = 0, it follows that both satisfy \z\ < 1 for all r E U . The case $ > 1 is 
similar. □ 

The following lemma is closely related to a similar result for the non-self- adjoint 
Anderson model in [TIT, Theorem 11]. 

Lemma 20 // c G fio- o-nd A G C then the space of all solutions of A^f = Xf is 
two-dimensional. If c is periodic with period p then the asymptotic behaviour as 
n — )■ ±00 of the solutions is determined by the solutions zi, Z2 of the polynomial 
z^ — r(A)2; + 7 = 0, where r(A) is a monic polynomial in A with degree p, given 
by r(A) = tr(Tp), where Tp = XpXp^i . . . Xi and 



X, 



n 



1 

-On A 



and 7 = det (Tp) = ztcr^. Ordering the two solutions so that \zi\ > \z2\, there 
are three cases: 
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1. 



A lies in the closed set 



Bc — {\:\zi\ = l and \z2\ = cr^}. 

This set is the spectrum of Ac, equivalently, the set of X for which Acf — 
Xf has a bounded solution. 

2. X lies in the open set 

Ic = {X:l> \zi\ > \Z2\ > aP}. 

This is the case if and only if all solutions of Acf — Xf decay exponen- 
tially as n ^ +00. 

3. X lies in the open set 

Oc = {A : l^^il > 1 and \z2\ < cr^}. 

This is the case if and only if there exists a solution of Acf — Xf that 
decays exponentially as n ^ +00 and grows exponentially as n ^ —00, 
and another solution that decays exponentially as n ^ —00 and grows 
exponentially as n ^ +00. 

Proof The sequence / : Z — > C is a solution of Acf — Xf if and only if 
fn+i + Cnfn-1 — Xfn for all n G Z. This recurrence relation can be rewritten 
in the form 




where T^ — . . . Xi. If c is periodic with period p, then the asymptotic 

behaviour of the two-dimensional space of eigenfunctions / is determined by 
the magnitude of the eigenvalues zi, Z2 of Tp. These are the solutions of the 
equation z"^ — tz-\-^ — Q where r = tr (Tp) and 7 = det(Tp). A simple induction 
estabhshes that the (i, j)-th entry of Tp is a polynomial in A with degree less 
than p unless z = j' = 2 in which case it is a monic polynomial with degree p. 
Therefore r is a monic polynomial in A with degree p. However 

p 

det(Tp) = JJ det(X^) = ci . . . Cp = 

r=l 
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does not depend on A. The continuous dependence of the roots of a polynomial 
on its coefficients implies that Be is closed while Ic and Oc are open. An 



application of Lemma 19 now completes the proof. One sees, in particular, 
that 

Spec(A,) = Ee = {A:$(r,7) = l}. 

□ 

Our next lemma enables us to determine the sets Ic and Oc for certain impor- 
tant periodic sequences c, and to determine the spectra of certain paired peri- 



odic operators. We continue with the assumptions and notation of Lemma 20 



Lemma 21 If V is a connected component of C\Bc then V Ic or V ^ 
Oc- If V is unbounded then V C Oc, and if E V then V C J^. // C\Bc 
has exactly two components then the bounded component equals Ic and the 
unbounded component equals Oc- 

Proof We first observe that V, Ic and Oc are all open sets and that their 
definitions imply directly that Ic, Oc are disjoint. Therefore V = {V (1 Ic) U 
(y n Oc), where the two intersections on the right-hand side are disjoint. Since 
V is connected, it follows that V = V O Ic or V = V H Oc This completes the 
proof of the first statement. 



Lemma 20 case 1 implies that 

B^ = Spec(A,) C {A : |A| < 1 + a}. 

Therefore C\Bc has only one unbounded component V and it contains {A : 
|A| > 1 + cr}. To prove that V ^ Oc it is sufficient by the first part of this 
proof to find a single point X E V (1 Oc- The fact that r is a polynomial with 
degree p implies that |t(A)| — )■ oo as |A| — )■ oo. This implies that the solutions 
of z"^ — t{X)z -|- 7 = 0, where 7 = icr^, are z ~ r(A) and z ~ 7/r(A) to leading 
order for all large enough |A|. Therefore A G Oc for all such A. 

The proof is completed by proving that G Ic- For A = one has Tp = 
XpXp-i . . . Xi where each Xr is of the form (^^^ J If p = 2m it follows that 

Tp = *Q ^°,„ j . The fundamental equation must therefore take one of the 
forms - 2(t"z + a^"^ = 0, + 2a"'z + a^'" = or - a^™ = 0. In each 
case both solutions have modulus cr^/^ < 1. The same holds ii p = 2m + 1. 

The final statement of the lemma follows from the following observations. 
There must be a component of C \ -Be that contains and there must be 
an unbounded component. The first part of the proof shows that these are 
distinct, and the extra hypothesis is that there are no other components. □ 

Our next task is to determine the sets Be, Ic and Oc for certain particular 
periodic sequences. 



16 



Lemma 22 If Cn = (J for all n & Z then Spec{Ac) is the ellipse 



Spec{Ac 



u + iv 



u 



+ 



" ^1 + ct2-2(tcos(2^) 



(27) 
(28) 



Moreover the interior U of the ellipse equals 1^ and the exterior V of the ellipse 
equals Oc- 



1 
-a X 



, SO r(A) = A and 7 = cr. Using (25) 



Proof We have p = 1 and Ti = 
we deduce that Spec(ylc) is given by (27). The proof is completed by using 
Lemma [2T1 □ 



Lemma 23 //c„ = —a for all ?2 G Z then Spec(Ac) is the ellipse 



Spec(y4c) 



u + iv : 



pe : p 



u 



+ 



V 



(29) 
(30) 



/I + (t2 + 2(7 cos(2^) 

Moreover 1^ is the interior of the ellipse and Oc is the exterior of the ellipse 



Proof We have p = 1 and Ti : 
rest of proof, which is almost identical to that of Lemma [22 



1 



, so r(A) = A and 7 = —a. We omit the 



□ 



In the following lemma, starting from Lemmas 22 and 23, and making succes- 
sive applications of Lemma 18, we compute the spectra of a family of periodic 
sequences, namely the sequences = ac^"'^^ G fio-, defined by (14). By 
Lemma 14, these sequences are periodic of period < 4". 

This next lemma applies for < cr < 1. The corresponding result for o" = 1 is 
equations (17) and (18) above. 

Lemma 24 Suppose n G Z+ and = ac^"''^\ c~ = ac^"'~\ Then 

Spec{Ac^) = {pe'' ■.p = p^{e,a)} (31) 

where 



Po ^) 



(l + a2-2(Tcos2^)^/^' 
and, for n G N, 



Po or J 



P^(^^,^)=(P^(2'^^,^'")) 



1/2" 



1 - O 



l + a2 + 2a cos 2^)^/^^ 

n + l\V2" 



(1 + T 2a2" COS (2"+i^)) 
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1/2" 



Moreover, 
and 



lc± = {pe'':0<p<ptie,a)} 
0,^ = {pe^':p>pi{e,a)}. 



Proof Our proof of (31) is by induction. We note first that (31) holds for 



n = by Lemmas 22 and 23 Suppose now that (31 ) holds for some n > and 



all < < 1. Then 



Spec(A^.,(„,±)) = {pe^' : p = p^i9,a')} = { p^' : p = {pt{2-9,a'-) 



1/2" 



Further, since ac^"^^^'"^^ = crr+(c*^"'+'') = rcr,+ (o"^c^"''^''), it follows from Lemma 
mthat 

Spec(A^c(„+i,±)) = |±V^ : A G Spec(y4^2c(n,±))| . 
Combining these equations, we see that (IMl) holds with n replaced by n + 1. 



Thus (31) follows by induction. 



The formulae for Ic± and Oc± follow from (31) and Lemma 21 



□ 



We remark that p„ {6, a) = p^{6 ±7t/2'"'^^, a), so that the spectra of Ac± in the 
above lemma are related by 

Spec(v4,+ ) = e±'"/2"+i Spec(A,-). 

This is a symmetry which is surprising from an inspection of the sequences 
c^, which need not even have the same period. (For example, as observed in 
Section |4| c+ has period 4 and c~ period 2 in the case n = 1.) 

In principle, since is periodic, (31 ) should be computable alternatively from 



the characterisation of the spectrum for general periodic sequences in Lemma 
20 As an example of this, for the sequence c~ = ac^^'~^ which has period 2, 
with c~ = (— l)"cT, the transfer matrix T2 is given by 



X2X1 



1 

-a A 



1 
a A 



a A 
o-A -a + A^ 



Applying Lemmas 19 and 20 with r = A and 7 = — a , we find that Spec(74c 



is the set of all A = m + iv for which 



- i;2)2 (2uvy 



2^2 



If one puts A = pe*^, then this may be rewritten in the form (|3l|). 
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The main point of the above theory and calculations are to prove and pre- 
pare the use of the following result on operators that are paired periodic 
operators. To state this result let us introduce the notations 

^»={- + '!'^^+(r^<l} (32) 

and 



so that and are the interiors of the ellipses introduced in Lemmas 22 
and I23I 

The following lemma is analogous to [TO^, Theorem 12], proved there for the 
non-self- adjoint Anderson model. 

Lemma 25 Suppose that c E ila is periodic and r G {a,— a}, and define 
c* G Qa by c* = Cn, for n > 0, and c*^ = r for n < 0. Then 

Spec(Ae*) DTc\Er. 
Proof Since Spec(Ac*) is closed it is enough to show that Spec(Ac*) D Ic\Er. 
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smce 



So suppose that A G Ic\Er. Then, by Lemmas 22 and 23 and Lemma 
A ^ Et-, it follows that there exists a non-trivial solution / of Ac*f = A/ such 
that /n — ^ exponentially as n — )■ — 00. Since \ E Ic, again applying Lemma 



20, it follows that this solution / also decays exponentially as n — )• +00. Thus 



A is an eigenvalue of Ac* so A G Spec(y4c.). □ 

7 First proof of the main theorem 



This section is devoted to the proof of Theorem 27, in which we assume that 
< a < 1. 

Lemma 26 We have 

Spec(AJ C |1 - a, 1 + al (34) 

and 

Spec{Ac) C{x + iy:\x\ + \y\< ^/2^^Ta^}, (35) 
for every choice of c EQa- 
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Proof We regard VcR as a small perturbation of L in the identity Ac = VcR+L, 
noted in the proof of Lemma |5} Since L is a unitary operator with spectrum 
{z : \z\ = 1}, we have 

ii(^-^/)~'ii = ii-kir' 



for all z not on the unit circle. The inclusion (34) now follows from ||K;-R|| = (J 
by a perturbation argument; see pT| Th. 9.2.13]. 

The inclusion (35) depends on an estimate of the numerical range of Ac. Fol- 
lowing [m Section 9.3], x + iy G Num(Ac) if there exists / G £^(Z) such that 
]|/]| =1 and X + iy = {Acf, /). This implies that 

X = liiAc + ADf, /), y = -'-{{Ac - ADf, /). 

Therefore 

x + y=^{BfJ) 

where 

B = iAc + A:)-^iAc-A:). 
A simple calculation shows that Bm,n = unless \m — n\ = 1, while 



Bn,n+1 = Bn+l,n = (1 ± (x) - i{l =F O") 



Therefore |-B„+i^„| = |-B„,„+i| = v2(l + cr^) for all G Z and 

x + y< ]^\\B\\ < ^2(1 + ^2). 

The other three steps in the proof of the bound for \x\ + \y\ are similar. □ 
The statement of our main theorem refers to the open set 

H^ = E^n E^„, (36) 

the intersection of the ellipses and E_^. This set satisfies 

|0,l-al c:^C |0,r.] (37) 

where 



2 



I- a 

T = — — 

Theorem 27 // < a < 1 then 

{A : |A| < 1}\H„ C S„. 



(3^ 
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Proof Note first that if and p^(6', a) are defined as in Lemma 24, then 



1-a 



\ 1 /2" 

2n + l \ 



l + (r 
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for all 6* G R, so that Ic± D {A : |A| < Pa,n}- Thus, defining c* G as in 



Lemma [25| with c = c"*" or c and r = ±cr, we see from Lemma |25j that 

Spec(A,.) dZ\E,D {X:\X\< p,,„} \ E,. (39) 
Applying Proposition [T], it follows that, for all n G N, 

S.D{X:\X\<p^^n}\H,. 
The theorem follows since sup„ Po-,n = 1 and 5*0- is closed. □ 





Figure 1: Plots of Spec(j4c) for the case when c is periodic and a = 0.5. The two 
plots show the sets ttat^o- C S^, the union of the spectra for ah sequences c of period 
< A^, for N = 2 (left) and = 12 (right). The two ellipses visible in the left-hand 



plot, the boundaries of Ef^ and E^a- defined in (32) and (33), are the components 
of TTi^o-; the other closed curve is Spec(Ac) for c = ac^^'~\ i.e. Cn = (— 1)"ct, given 



explicitly in Lemma 24 Also shown in each plot are the boundaries of the inclusion 
namely the circles of radius 1 it a, and, in dashed lines, the 



sets from Lemma 26 



boundary of the set (35). The boundary of H^, which we conjecture is a hole in the 
spectrum 5o- , is high : 



ighted in a thicker line. 



Lemma 26 and Theorem |27] together establish that there is a hole in S^^ which 
is at least as big as {A : |A| < 1 — cr} and which is no larger than H^. The 



21 



numerical computations we have been able to carry out are consistent with a 
hypothesis that the hole is precisely the set i.e. they are consistent with 
a hypothesis that Spec (Ac) fl H^j = for every c G and hence for every 
c e £a. 

It should be pointed out, however, that these numerical computations are only 
for instances where c is periodic, for which we have a characterisation of the 
spectrum in Lemma 20 Thus, strictly speaking, our calculations are evidence 
of the possibly weaker result that tTqcct H H^j = 0; they become evidence that 
Spec(ylc) n i^o- = with a hypothesis that TCoo,a is dense in the part of S'o- that 
is contained in the unit disc. This latter statement may or may not be true for 
a G (0, 1), but we have shown in Theorem 15 that it is true for a = 1. 

As an example of the numerical computations we have carried out, the right 
hand side of Figure [T] shows the union of Spec(Ac) over all periodic c G l^cr for 
which the period < 12. It is clear from this figure that 7ri2,o- H Her = for 



a = 0.5. We note that, rather than using the characterisation in Lemma 20 



we use for these computations the standard Bloch-decomposition formula (e.g. 
tllj) that 

Spec(A,) = U Spec(4JP<'^)), (40) 
l"l=i 

where Ai'^'^'^'^^ is the N x N matrix defined in (51) below. 




Figure 2: Plots of 10^ instances of Spec(74c,o'''™ ) for a = 0.5, computed as described 
in the text, with N randomly chosen in the range 1 < < 100. 
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It is not feasible to calculate 'Kn,cti the union of all 2 periodic spectra of 
period A^, for very much larger values of A^. In Figure |2] we sample 7rioo,cr; 
for a = 0.5^ plotting the union of the spectra of 10^ randomly chosen N x N 
matrices aJ^'^'^^\ By randomly chosen we mean here that, in each realisation, 
N e {1, 100} is randomly chosen, with higher probabilities for the smaller 
matrix sizes, and then the vector c = (ci, ca?) is randomly chosen, with each 
Cn = ±cr independent and identically distributed with Pr(c„ = cr) = 0.5, and 
finally the phase factor a is randomly chosen from a uniform measure on the 
unit circle. We see in the figure that clearly, as they have to, the spectra are 
constrained to lie in the inclusion sets shown in Lemma |26l We also note that 
all the spectra lie outside Hu. 



8 Second proof of the main theorem 



In Theorem 34 of this section we provide a second proof that 

{A : |A| < 1}\H, C (41) 

for all cr e (0, 1). This proof has a lot in common with the previous one, but 
it reveals more about the asymptotic behaviour of the solutions of the second 
order recurrence relation for certain choices of c G fia- A key role in this section 
is played by the sequence Ce G fl^ defined in Section |4| which sequence was 
central to the proof of Theorem [T] that appears in [5J (see Theorem 12 above 
for more details). 

We start with some calculations that do not depend on a. Throughout this 
section Cn G {±1} is defined for all n > 1 by the rules Ci = 1, C2„ = C2„_iC„ 
and C2„ + C2n+i = 0. (In other words, Cn = Ce^n, for n > 1.) The first few 
values are shown in Table [ij We will obtain a bound on a transfer matrix Tm^x 



associated with this sequence and use this bound to prove Theorem [34 

Let M : Z+ — 7- C be the solution of Un+i = Xun — CnUn-i such that uq = 
and Ml = 1, so that u is the restriction to Z4. of the bi-infinite sequence Ue 
already studied in |5j and discussed in Section |4j We have observed already 
in Proposition [TT] that u„ is a polynomial of degree n — 1 in A with integer 
coefficients for all n > 2. Similarly, if v : Z+ — )■ C is the solution of Vn+i = 
Xvn—CnVn-i such that f = 1 and vi = 0, it is easy to see that f„ is a polynomial 
of degree n — 2 with integer coefficients for all n > 2. 

One may check the computations of Um and Vm in Table [T] by confirming the 
determinantal identity UmVm+i — VmUm+i = ±1 for all m > 1, the left hand 
side being a polynomial in A. Here we are referring to a determinantal identity 
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n 






Un 




Vn 


1 


1 




1 







2 


1 




A 




-1 


3 


-1 




- 1 




-A 


4 


-1 




A3 


-A2 


- 1 


5 


1 




- 1 


-A3- 


-2A 


6 


-1 


A^ 


-A 


-A^ - A^ 


+ 1 


7 


1 


A^ + A^ 


- 1 


-A^ - 2A3 


-A 


8 


-1 




A^ 


-A^ - A^ 


- 1 


9 


1 


A^ + A^ + A^ 


- 1 


-A^ - 2\' - 2A3 - 


-2A 



Table 1: Values of c„, f „ for 1 < n < 9. 



for the transfer matrix Tm,\, defined by 



which transfers the data of any solution of Xn+i = Xxn — CnXn-i from {0, 1} to 
{m, m + 1} in the sense that 



It is easy to verify (see Lemma 20 above) that 

Tm,\ = X„iXrn-l . . . Xi (42) 

where 

^' = (-11) («' 

has determinant cv G {±1} for every r > 1. 

The proof of Theorem [34] we will give below was motivated by numerical evi- 



dence that 

tr(T2.,A) = A^'- - 2 

holds for all r > 1 and all A G C; see Table [2j We prove this crucial identity 
in Lemma mi 

Lemma 28 Let T be a 2x2 matrix with determinant 6 and trace r. Ifr'^^AS 
and 7 is the absolute value of the larger root of — tz + 5 = then there exists 
a constant b such that 

IIT^II < bY 
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n 




) - 


Un+l + Vn 


1 






X 


2 






A2-2 


3 






A3 -A 


4 






A^-2 


5 




X' 


- A3 - 3A 


6 






A^-A^ 


7 






- 2A3 - A 


8 






A«-2 



Table 2: Values of tr(T„,A) for 1 < n < 8. 



for all r > 1. If r'^ = 46 then for every e > there exists a constant 6^ such 
that 

\\T'-\\<Ki^ + EY (44) 

for all r > 1. 

Proof The eigenvalues z± of T are the roots z of z^ — tz + 6 = 0. The condition 
7^ 45 implies that the eigenvalues are distinct, so T is diagonalizable - there 
exists an invertible matrix B such that 

T = Bi ^+ ^ ]b-\ 



Therefore 



IT'II = \\B ( 'q+ ) 5-^11 < \\B\\ ||5-^||7^ 



The slightly worse bound (44) is obtained when = 46 because one has to 
use the Jordan canonical form for T. □ 

Lemma 29 The identity det(T2n = 1 holds for all n > 1 and all X E C. 

Proof If m G N then ^ and disl) imply 



det(T2„,A) = lldet{X2rX: 



2t-1, 

r=l 
m 



W{c2rC2r-l] 

r=l 
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r=l 



det{Tm,x)- 



It follows by induction that 

det(T2n_A) 



det(ri,; 



Cl 



1. 



□ 



The following lemma depends on Proposition 11 above, abstracted from [5], 
which notes properties of the integer coefficients pij of the polynomials 



Ui 



j=i 



Lemma 30 The polynomial Um is even for odd m and odd for even m. Its 
leading term is A'""^. //m = 2" and n > 2 then 



Ur 



m/4 

_l + A"^/2^a,A 

r=0 



2r 



(45) 
(46) 



where aj. G {0, 1, —1} for all r. 

Proof The statements in the first two sentences may be proved by induction, 



using the definition of Um- We prove (45) and ( |46| ) for m = 2" by induction in 
n, noting that both hold for n < 3; see Table [T] As in Proposition 11, let Y 
denote the set of (z, j) G such that pij ^ 0. 

To prove (45) suppose that (2"+^, j) G Y . Proposition [TT| implies that j is even 
and that (2",j/2) G Y . The inductive hypothesis now implies that j/2 = 2", 
so j = 2"+^ 



To prove (46) suppose that (2""*"^ + l,j) G Y . Proposition 11 implies that j is 
odd and either (2^^, (j + l)/2) G F or (2" + 1, (j + l)/2) G Y. In the first case 
we have already proved that (j + l)/2 = 2", so j = T^^^ — 1. In the second 
case the inductive hypothesis implies that (j ' + l)/2 > 2"~^ + 1, so j > 2" + 1 
or (j + l)/2 = 1, so J = 1. 

We finally have to evaluate the constant coefficient 7^ of Um when m = 2" + 1. 
This may be done by considering the defining recurrence relation in the case 
A = 0, namely 7^+1 = — Cr7,._i subject to 70 = and 71 = 1. This implies that 



72m+l 



m 
r=l 
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for all m G N. Therefore 75 = 73 = — 1 and 

m 

78m+l = W (c8rC8r-2C8r-4C8r-6) 



r=l 

m 



— W (c4rC8r-lC8r-2C4r-2C8r-5C8r-6) 
r=l 
m 

= W (c4rC4r-2) 

r=l 
= 74m+l 

for all m G N. A simple induction now implies that 7^ = — 1 for m = 2" + 1 
and all n > 1. □ 

Lemma 31 If m = 2" and n >2 then 

T = tr(T„,A) =Vm + Um+l = A"* - 2 (47) 

for all X e C. 



Proof The proof uses the identity VmUm+i — UmVm+i = 1 of Lemma 29 to- 
gether with the two identities proved in Lemma 30 These are identities within 
the commutative ring Z(A) of all polynomials with integer coefficents in the 
indeterminate quantity A, but they imply similar identities in the commuta- 
tive ring Z(A) of all polynomials with integer coefficients in an indeterminate 
quantity A that satisfies the identity A™""^ = 0. (Equivalently one may start by 
disregarding all terms in the identities that involve A^ with r > m — 1.) The 
identities then simplify to 

where 

m/4 

p(A) = A™/2^a,A2r 

r=0 

satisfies = in Z(A). The second equation can be solved for Vm, yielding 

Vm = -1-p 

and hence r = — 2. Returning to the original variable A one deduces that 

T = -2+ J2 P^y. 

r>m— 1 



But (cf. Lemma 30 ) it is easily shown that f „ is an even polynomial of degree 
n — 2 for odd n. Thus, and by Lemma 30 , r = Vm + Um+i is an even polynomial 
of degree m with leading coefficient 1, so r = A™" — 2. □ 
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Lemma 32 Following the assumptions and notation of Lemma 28_, suppose 
that 5 = 1 and that there exist m G Z_|_ and /i G C such that t = fi"^ — 2. Then 
there exists a constant b such that 



IT'W < 64'^max(|/i|'^"^,l) 



for all r > 1. 



Proof 

Case 1. If < 1 it suffices to obtain bounds on the solutions z± of z"^ — tz + 
1 = when |r| < 3. The solutions satisfy 



± 



< 



< 4. 



Lemma 28 now implies that ||T''|| < 64'" for all r > 1. 



Case 2. If > 1 it suffices to obtain bounds on the solutions z± of z"^ — rz + 
1 = when |r| < Sl/xl'". The solutions satisfy 



± 



< 4|/i| 



Lemma 28 now implies that IIT*"!! < 64'"|/i|'"™ for all r > 1. 



□ 



Lemma 33 Let X„, n & Z, be a periodic sequence of 2x2 matrices with period 
m and let Tr = X^Xr-i . . .Xi for all r > 1. If there exist constants bo, 7 such 
that ||(Tm)*|| < boY f(^f^ 0.^^ s > 0, then there exists a constant 62 such that 
\\Tr\\ < 627'^/™ for allr>l. 

Proof Every r G Z_|_ may be written in the form r = sm + v where s > and 
< f < m. Using the identity = {TmY, one obtains 

||X^|| ||X.f.X,, — I . . . ^sm+l-^smll 
= ||Xt,Xy_i . . . Xi(Tm)'^|| 

< ||X,X,_i...Xi||6o7'™ 

< bob^Y 

< 627'/'", 

where 62 = ^o^i and 

61= max {||X,X,_i...Xi||}. 

0<D<m— 1 



□ 
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Theorem 34 One has 

{A : |A| < 1}\H, C (48) 

for all a G (0, 1). 

Proof Given a G (0, 1) we put m = 2'^ where d G N is large enough to yield 

<h = 4-^/". (49) 

We use the identities 



5 = det(T, 



m,fj,J 



and 



r = tr(T™,^) = /i'- - 2 

proved in Lemmas [29] and [31] and valid for all yU G C. Let c G i^^r be the 
periodic sequence with period m such that = aCn for all 1 < n < m. The 
main task is to prove that if |A| < h then all solutions ^ : Z — > C of 



^71+1 — A.^„ — Cn^n-l 



(50) 



decay exponentially as n — )■ +oo. This will imply, by Lemma [20} and using 
the notations of that lemma, that 

J, D {A : |A| < h}. 



Arguing as in the proof of Theorem |27[ it will then follow from Lemma 25 and 
Proposition [T] that 

S.D{X:\X\<h}\H,, 
this holding for any h = 4"^/"* such that (49) holds and m = 2'^, so that 

D {X : \X\ < 1} \ H^. 



Since 5*0- is closed, (48) will follow. 



Thus it remains only to show that all solutions of ( 50 ) decay exponentially at 
+00. To see that this holds, define x„ = c7~"'/^^„ and = a~^^'^X so that (50) 
may be rewritten in the form 



for 1 < n < m. Where 9 = max(l, |/i|). Lemma 32 now yields 



< 64"^^^ 



for all r G N. Lemma |33| with 7 = 46'™ implies 
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and hence 

\xr\ < h4r^"'9^, 
again for all r e N. Hence, where = max((T^/^, |A|), 

for all r e N. Since < < /i, it follows that ^ decays exponentially. 



□ 



9 Semi-infinite and finite matrices 



All our results so far have focused on calculations of the spectrum of the bi- 
infinite matrix Ac. In this final section we say something about the spectrum 
of the semi-infinite matrix 



A1: 



( 1 

ci 



1 

C2 



V 



/ 



in the case that c = (ci, C2, ...) G {icr}^ is pseudo-ergodic (contains every finite 
sequence of icr's as a consecutive sequence). We also say something (though 
have mainly unanswered questions) about the finite N x N matrices 



/ 1 

ci 1 

C2 

V 



cn-1 / 



and 



/ 1 

ci 



1 

C2 



(51) 



cn-1 / 



Here Ac^^ is tridiagonal, A^^''^'^^^ is tridiagonal except for "periodising" entries 
in row 1 column N and row N column 1 (in these entries we assume that 
\a\ — 1), and each Cj — icr: we have in mind particularly the random case 
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where the Cj's are independent and identically distributed random variables 
taking the values ±a. 

Our main result on the spectrum of A^., proved in the previous sections, is 
that it contains the set {A : |A| < 1} \ H„. We suspect that H„ is a genuine 
hole in the spectrum for < cr < 1, i.e. that H„ fl 6*0- = 0. We have not 



shown this result but have shown in Lemma 26 the weaker result that {A : 
|A| < 1 — a} n Sa = 0. Our first result in this section is that this hole is 
not present in the spectrum of the semi-infinite matrix. The proof depends 
on recent results on semi-infinite pseudo-ergodic operators due to Lindner and 
Roch [TU] , derived using characterisations of the index of Fredholm operators, 
whose matrix representations are banded semi-infinite matrices, in terms of so- 
called "plus indices" of hmit operators, these characterisations derived using 
i^-theory results for C*-algebras in [20] . 



Theorem 35 Suppose c G {±cr}^ is pseudo-ergodic. If a = 1 thenSpec{A^) = 
S. For all a E (0, 1], {A : |A| < 1} C Spec(A+). 

Proof In the case that a = 1 it is shown in [6] that Spec(y4^) = S. Thus, for 

(7=1, 

{A : |A| < 1} C 5 = Spec(A+) 

follows from Theorem [T] (or [5l Theorem 2.3]). For all a G (0, 1] it follows from 
[T9| Theorem 2.1] that the essential spectrum of A^, i.e. the set of A G C for 
which A'^ — A/+ is not Fredholm (here is the identity operator on £^(N)), 



is the set So-- Thus and by Theorem 27 



({A:|A|<l}\i7,)cS. cSpec(A+). 

It remains to show that ifo- C Spec(A+). But, applying [IQ,, Theorem 2.4] 
(note that the set E-{U,W) in the notation of [T9l Theorem 2.4] is precisely 
the set H„ for this operator), it follows that, for A G i/o-, either — A/+ is 
not Fredholm or A^ — XI^ is Fredholm with index 1: in either of these cases 
A G Spec(A+). □ 

Our other result in this section is to say something about the spectra (sets 
of eigenvalues) of the finite matrices Ai^^ and A^^'^^'^K The notations tttv.o- 
and defined in and above equation ([T]) (and ttn and tToo are our 

abbreviations for a = 1). 

Theorem 36 // < (x < 1, |a| = 1 and c G {±(t}^ , then 

Spec(/l2'P^'"^) C vr^v,. C tt^,^ C 

while 

Spec(A^^^) C \/a7T2N+2 C y/anoo C y/aS. 
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If X = X + iy is an eigenvalue of A[^'^'^^^ then 1 — cr < |A| < 1 + cr and 



\x 



+ \y\ < v^2(1 + (t2), while ifX is an eigenvalue of Ai^^ then \x\ + \y\ < 2a/(t. 



Proof The first of these statements is clear from the definition of TiN,a, (40), 
and Proposition [l] which gives that 7ioo,a C S^- The second of these statements 
is shown for a = 1 in [6, Theorem 4.1]. The second statement follows for 
< (T < 1 by the observation that, where d G {±1}^, c = ad E {±0"}^, 
and is the diagonal matrix with leading diagonal (1, cr^/^, a, cr*^^~^^/^), it 
holds that D^A^i^^D^ = y/a A^J^\ The last sentence then follows from Lemma 

□ 



Note that in the last sentence of the above theorem the condition + < 2y/a 
implies both that |A| < 1 + a and that |x| + \y\ < \/2{l + a^) . 




Figure 3: Plots of Spec(^^^^) (left) and Spec(^SP'"'^) (right) for a case when 
N = 2000, a = 0.9025, a = 1, and the entries of the vector c = {ci, ...,C]\f) are 
independent and identically distributed with Pr(cj = ±a) = 0.5 for each j (the same 
vector c is used in the two plots). 

In Figure [i] we plot the spectra of Ai^'^ and Ai^'^^^'^ for = 2000 and a = 1 
for a typical realisation with the entries c G {±a}^ randomly chosen with the 
Cj independently and identically distributed with Pr(cj = a) = 0.5 and a = 
0.9025 so that ^/o' = 0.95 (the several other realisations we have computed are 
very close in appearance to these plots). Theorem 36 tells us that Spec(Ai^'') C 



0.957roo C 0.955* and that Spec(74c^''"^'^-') C 5*0.90255 and that if A = a; + iy is an 
eigenvalue of Ai^^ then |x| + \y\ < 1.9, while if A is an eigenvalue of Ai^"^"^"^"^ 
then 0.075 < |A| < 1.9025 and |x| + \y\ < ^^2(1 + ct2) ~ 1.905. 
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It is clear from Figure [3] that Theorem 36 is only the beginning of the story. 
We observe in the figure a hole in the spectrum of Ai^'^'^'^\ but it is a hole of 
radius approximately 0.6, not 0.075, with a large proportion of the eigenvalues 
positioned on the boundary of this hole, while outside the hole the spectra of 
aZ""'^ and ^ appear near identical. The same qualitative behaviour is vis- 
ible in Figure |4| which is a similar plot except that a is reduced to 0.5 and we 
change the probability distribution, making it twice as likely that the entries 
of the vector c are —a rather than a. This change of probability distribution 
introduces an asymmetry, in particular an asymmetry in the hole in the spec- 
trum (if we instead compute with Pr(cj = a) = 1/2 then typical realisations 
have spectra which are approximately invariant under the dihedral symmetry 
group D2 of the square). Of course our methods, which are not probabilistic, 
have nothing to say about such asymmetries, indeed nothing, beyond Theorem 
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to say about the almost sure spectra of Ai^"^ or Ai^^'^^^ as N ^ 00. 




Figure 4: Plots of Spec(ArO (left) and Spec(^^;a''''"'0 (right) for a case when 
N = 2000, C7 = 0.5 and the entries of the vector c = (ci, cat) are independent and 
identically distributed with Fv{cj = a) = 1/3. 
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